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FOREWORD 

The  current  work  on  variable  moduli  models  is  part  of  Contract  No. 
DACA-39-67-C-00li8,  "Investigation  of  Ground  Shock  Effects  in  Nonlinear 
Rysteretic  Media,"  being  conducted  for  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  under  DASA  sponsorship. 

The  contract  is  being  monitored  by  Dr.  J.  S.  Zelasko,  Impulse  Loads 
Section,  Soil  Dynamics  1  ranch,  under  the  general  supervision  of  Mr.  A.  A. 
Maxwell,  Acting  Chief,  Coils  r'ivision,  WES .  COL  Levi  A.  Brown,  CE,  is 
the  contracting  officer. 


ABSTRACT 


The  problem  of  one  dimensional  plane  wave  propagation 
is  considered  in  isotropic,  nonlinear,  hysteretic  materials. 
Solutions  for  surface  loads  of  step  pressure  and/or  shear, 
in  cases  of  loading  or  unloading  are  studied,  and  the  presence 
of  shocks  and  regions  of  continuous  stress  change  are  discussed. 

The  mathematical  models  used  are  a  generalization  of 
conventional  elastic  models.  The  moduli  K  are  assumed  to  be 
functions  of  tht  stress  invariant  J  ^  ,  and  the  moduli  G  are 
functions  of  both  and  the  stress  deviators.  Different 
expressions  for  K  and  G  are  used  during  loading  and  unloading, 
leading  to  energy  loss  through  hysteresis. 
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LIST  OF  SYMBOLS*^ 

Shear  modulus,  shear  moduli  in  loading 
and  unload ing , 

Elastic  shear  modulus,  elastic  bulk  modulus. 
Nondimensional  shear  modulus. 

Coefficients  used  in  variable  shear  moduli, 
Eqs.  (24)  and  (25). 

Stress  invariants. 

Bulk  modulus,  bulk  moduli  in  loading  and 
unloading. 

Nondimensional  bulk  modulus. 

Coefficients  used  in  variable  bulk  moduli, 
Eqs.  (22)  and  (23). 

Fractional  coefficient. 

Stress  deviators. 

Time . 

Nondimensional  dependent  variable  defined 
by  Eq.  (14). 

Velocity  components  in  x,  y,  z  direction, 
respectively. 

Cartesian  coordinates. 

Volumetric  strain. 

Deviatoric  strains. 

Mass  density. 

Normal  stress  components. 

Shear  stress. 


Differentiation  with  respect  to  time. 
Differentiation  with  respect  to  U. 


Other  symbols  are  defined 


as  they  appear  in  the  text. 
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I  INTRODUCTION . 


The  following  is  a  study  of  one  dimensional  wave  propagation 
of  pressure  and/or  shear  for  isotropic  materials  with  nonlinear 
properties  resulting  in  hysteretic  effects.  The  one  dimensional 
plane  case  is  considered  because  it  permits  easy  mathematical 
treatment  leading  to  an  understanding  of  the  behavior  in  multi¬ 
dimensional  situations.  In  addition,  it  can  be  used  as  a  check 
on  numerical  schemes  and  gives  limiting  solutions  for  two 
dimensional  situations  at  high  Mach  numbers.  Solutions  for  the 
one  dimensional  problems  for  the  step  loads  studied  exist  for 
all  input  combinations  of  pressure  and  shear,  and  are  unique, 
giving  at  least  a  hint  that  the  artificially  constructed  des¬ 
cription  of  the  material  does  not  lead  to  absurdities. 


The  mathematical  models  used  are  a  generalization  of  con¬ 
ventional  elastic  models,  the  elastic  moduli  K  being  now 
functions  of  the  first  stress  invariant  J^,  ,  while  the  moduli  G 

are  functions  of  the  stress  invariants  and  J 2  .  The  functions 
K  and  G  differ  appropriately  during  loading  and  unloading, 
leading  to  energy  loss  through  hysteresis.  Following  Ref.  [1] 
the  governing  volumetric  relations  for  initial  loading  are 


J1  '  3KLD<J1>  Skk 


while  for  unloading  and  reloading 


J1  =  3KUN^J1)  ^kk 


(1) 


(2) 
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where 


K  >  K  ■>  0 
UN  -  LD 


Ic  is  noted  that  Eq .  (1)  permits  determination  of  as  a 
function  of  e  »  so  that  the  equation  can  be  written  in  the 
alternative  form  without  change  of  meaning. 

The  alternative  form  is  used  for  computational  convenience  in 
Ref.  [11  because  it  results  in  simple  expressions  for  the 
velocity  of  compressive  shocks. 

The  deviatoric  relations  are,  for  loading  and  reloading*), 


Sij  =  2GLD(J3  ’  J2)eij 


and  for  unloading 


Sij  =  2GUN(J1  ’  J  ? ^ e i  i 


where 


G  >  G  >0 
UN  -  LD  - 


2  iJ 


The  choice  of  shear  modulus  is  based  on  the  sign  of  > 

Gld  being  used  for  J2  >  0,  while  G^.  is  used  for  j,,  <_  0. 

Note  that  this  paper  uses  tensile  stresses  and  strains 
as  positive  while  Ref.  [1]  uses  compressive  stresses  and 
strains  as  positive. 


Whether  GLD  or  Gy^  is  appropriate  to  use  for  the  modulus 

of  reloading  is  a  question  still  undergoing  discussion 

and  investigation.  This  pe.per  assumes  tentatively  that 

G  is  appropriate. 

L  D 
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The  equations  of  motion  are 


Dt 

3x 


1 

3 


(12) 

(13) 


The  five  differential  equations  (7)  and  (10)  to  (13)  in  the 

five  unknowns  s^  ,  ,  T,  u  and  v  define  the  problem.  The 

two  equations,  (8)  and  (9),  serve  solely  for  the  determi- 

s 

nation  of  s  and  s  ,s  =  s  =  -  ,  a  trivial  relation 

y  z  ’  y  z  2 

which  will  not  be  carried  along. 


Following  a  procedure  used  in  a  previous  investigation 
of  plane  waves  in  elas to-plas tic  materials,  Ref.  [2], 
dimensional  considerations  require  that  all  stresses  and 
velocities  depend  on  a  nondimensional  combination  of  para¬ 
meters  and  of  the  variables  x  and  t, 


U 


7 


(14) 


Introduction  of  the  variable  U  is  conveniently  combined  with 
the  definition  of  the  nondimensional  quantities 


o 

G  ( J 


Sx  =  G 


8  = 


i  -  V 


T  = 


K(J1) 


(15) 


(16) 


i 
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Introducing  Eqs.  (14)  to  (16)  into  Eqs.  (7)  and  (10)  to 
(13)  yields,  after  elimination  of  the  velocities  u  and  v, 

(U2  -  s)  T*  *  0 


A&  j'  _  s'  „  o 
9k  Ji  Sx  0 


1  u  1  * 

<3  -  Tk>  J!  +  S::  “  ° 


(17) 

(18) 
(19) 


where  primes  indicate  differentiation  with  respect  to  U. 

These  three  differential  equations  apply  in  regions  of  con¬ 
tinuous  solution,  i.e.,  provided  the  derivatives  are  finite. 

The  differential  equations  are  linear  with  respect  to  the 

I 

first  derivatives  J ^  ,  etc.,  and  homogeneous.  They  thus 

i  i  i 

permit  only  the  solution  J  ^  -  s^  =  T  =0,  with  ^  =  constant, 

s  =  constant,  T  =  constant,  unless  the  determinant  of  the 
x 

coefficients  of  Eqs.  (17)  to  (19)  vanishes.  Noting  that 
Eq .  (17)  depends  on  T,  but  not  on  the  two  other  variables, 
while  Eqs.  (18)  and  (19)  do  not  contain  T,  the  respective 
subdeterminants  may  be  considered  separately. 

Equation  (17)  leads  to  the  conclusion  that  T  must  be 
constant  in  any  region,  unless 

U2  -  g  =  0  (20) 

Equations  (18)  and  (19)  indicate  that  s^  and  will 
be  constants  unless,  in  the  region  considered,  the  determinant 
of  the  coefficients  vanishes. 

3U 2  -  3k  -  4g  =  0  (21) 


While  it  will  be  seen  in  the  applications  that  this  does 
not  actually  occur,  the  possibility  of  simultaneous  variation 
of  Eqs.  (20)  and  (21)  must  be  allowed  for. 

The  continuous  solutions  of  the  differential  equations 
(17)  to  (19)  do  not  completely  describe  the  situations  which 
may  be  en'ountered.  There  may  be  locations  of  discontinuities 
(shocks)  which  must  be  obtained  from  other  considerations.  It 
is  known  that  materials  with  a  hardening  pressure-volume 
relation  permit  (pressure)  shocks.  For  the  materials  considered 
in  this  study  the  velocity  of  such  shocks  may  be  obtained  by 
integrating  the  stress-strain  relations  in  uniaxial  strain. 

In  principle,  shear  shocks  are  also  possible,  but  such  shocks 
do  not  occur  for  the  dependency  of  G  as  a  function  of  con¬ 
sidered  here,  because  the  material  in  this  respect  is  a 
"softening  one". 

The  next  section  will  consider  specific  situations  for 
the  material  designated  Type  II  Variable  Modulus  Material, 

Ref.  (3].  The  pertinent  point  lies  not  in  the  details  of  the 
solution,  but  in  the  fact  that  there  is  always  a  solution, 
and  that  there  is  just  one  solution  even  if  the  possibility 


of  shocks  is  included. 


--  7  -- 

HI  SOLUTIONS  FOR  SPECIAL  CASES. 

Typical  suitable  expressions  for  the  moduli  of  soils  have 
been  obtained  for  uniaxial  and  triaxial  test  results  of  the 
type  represented  by  Figs.  2  and  3,  respectively.  The  ex¬ 
pressions  proposed  are  of  the  following  form. 


Case  A 


Shear  waves,  initial  loading. 


At  the  time  t  =  0,  a  stressless  half-space,  Fig.  1,  is 
subjected  to  a  load  Tq  ,  while  0^0  -  0.  Noting  that  Eq  .  (14) 
gives  U  =  0  when  x  =  0,  the  boundary  condition  on  the  surface 
is 


T(0)  =  To  (28) 

Further,  for  large  values  of  x  -*■  °°,  U  °°, 

T(“)  =  0  (29) 


If  a  continuous  solution  exists,  Eq.  (20), 


U 


2 


SLD 


(30) 


must  hold.  Noting  that  =  0,  s^  =  0,  Eqs.  (24)  and  (26) 
give 

U2  =  1  -  g x | T |  (31) 

This  simple  relation  defines  the  variation  of  t  as  a 
function  of  U  in  any  region  where  t  is  not  a  constant.  How¬ 
ever,  regions  where  T  is  simply  constant  must  also  be  in¬ 
cluded  in  the  construction  of  a  solution.  It  is  important 
that  Eq.  (31)  indicates  that  if  |t|  is  not  a  constant  it  will 
increase  towards  the  surface  where  U  is  smaller.  Equation 
(31)  is  therefore  suitable  for  the  particular  case  of  loading. 
Provided  g  ^  j  T  |  <  1  there  will  be  a  point,  =  it  where  the 
shear  stress  begins  to  increase,  see  Fig.  4a.  Equation  (3!.) 


__  9  -- 


describes  the  solution  until  the  point  U ^  ,  defined  by 


(32) 


is  reached.  For  U  <  c*ie  stress  T  simply  remains  constant, 
so  that  the  solution  shown  typically  in  Fig.  4a  is  obtained. 


The  relation  Eq .  (31)  being  monotonic,  it  is  impossible 
to  find  an  alternate  continuous  solution  to  satisfy  the 
boundary  condition*^  . 


No  solutions  containing  a  discontinuity  can  be  obtained 
because  the  weakening  character  of  the  stress-strain  relation¬ 
ship  in  shear.  Fig.  3,  does  not  permit  discontinuities  in 
shear . 


Case  B  -  Shear  waves,  unloading 


Let  the  half-space  be  initially  uniformly  stressed  by  a 

shear  stress  T  =  T  .  At  the  instant  t  =  0  the  surface  stress 

o 

is  removed.  The  boundary  conditions  are,  in  this  case, 


T (0 )  =  0 

T(w)  =  T  o 


(33) 


To  investigate  continuous  solutions  the  reasoning  of  Case  A 
applies  again,  but  the  expression  g^  of  Eq.  (25)  must  now 
be  used,  so  that 

°2  =  S7  +  g 4  |  T  ]  (34) 

'  It  is  also  not  possible  to  construct  alternative,  in¬ 
appropriate  solutions  of  the  type  discussed  in  Case  F. 
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and  gy  being  positive.  The  relation  indicates  that  jt|  in¬ 
creases  with  U,  as  required  for  an  unloading  case.  Figure  4b 

shows  the  solution  to  be  T  =  T  for  values  U  >  U,  ,  where 

o  1 

U1  -  S7  +  84To  (35) 

The  value  T  =  0  is  reached  for  =  /g y  and  T  remains  at  this 
value  for  U  <  • 

*  ) 

No  alternative  solution  involving  a  shear  shock  is 
possible,  because  of  the  "weakening"  stress -strain  relation 
for  unloading  in  shear  shown  in  Fig.  3. 

If  a  material  has  a  stress-strain  relation  for  unloading 
in  which  the  curvature  shown  in  Fig.  3  is  reversed,  a  shear 
shock  is  possible  and  would  give  a  solution  to  the  problem. 
However,  in  this  case  the  sign  of  g^  in  Eq.  (25)  would  be 
negative  and  there  would  be  no  continuous  solution  possible. 

The  uniqueness  and  existence  requirement  would  still  be 
satisfied. 

Case  C  -  Pressure  waves,  unloading . 

Consider  the  half-space  subject  to  a  uniform  state  of 
stress  defined  by  the  given  value  J.  <0  and  s  <  0,  so 

X  O  XO 

that  the  normal  surface  stress  is  a  =  s  +  4  J ,  <0.  At 

xo  xo  3  lo 

the  time  t  =  0  the  surface  load  is  reduced  to  a  fraction,  say 
n0xo  '  Tlie  boundary  conditions  in  terms  of  the  independent 


*) 


It  is  also  not  possible  to  construct  alternative,  in¬ 
appropriate  solutions  of  the  type  discussed  in  Case  F. 


variable  U  are 


sx(°) 

+  i  Jx(0)  =  nOxo 

(36) 

s  (°°) 

B  s  <0 

(37) 

X 

xo 

Jl<“) 

«  J,  <0 
lo 

(38) 

The  determination  of  continuous  solutions  is  more  com¬ 


plicated  than  in  Cases  A  and  B,  because  the  problem  concerns 

now  two  unknown  functions  J,  and  s  ,  while  the  earlier  cases 

1  x 

concerned  only  one  unknown.  As  previously  stated  the  differ¬ 
ential  equations  (18),  (19)  permit  no  other  solution  than 

J,  =  constant,  s  =  constant  and  thus  a  =  constant,  unless 
1  x  x 

the  determinantal  equation  (21)  is  satisfied.  If  this  is 
the  case,  Eqs.  (18)  and  (19)  are  inherently  equivalent  and 
only  one  of  them  need  be  retained.  It  is  convenient  here 
to  select  Eq .  (18) 


(39) 


which  relation  is  to  be  solved  in  conjunction  with  Eq.  (21). 


Before  considering  specific  cases  some  general  con¬ 
clusions  can  be  drawn  from  Eq.  (39).  The  quantities  g  and  k 

i  i 

are  inherently  positive,  so  that  and  s^  must  necessarily 
have  the  same  sign.  In  the  problem  of  unloading  considered 

I 

here  it  is  ttierefore  inherent  that  negative  values  of  and 
» 

s^  will  represent  unloading  as  long  «s  s^  <  0,  so  that  the 

unloading  relations,  Eqs.  (23)  and  (25),  for  k  and  g  are 
*) 

appropriate 


If  n  is  sufficiently  small,  s  will  change  in  sign  and 
Eq .  (24)  v ' 1 1  be  appropriate, x  For  simplicity,  th*s 

case  is  not  pursued  here. 
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In  order  to  find  out  if  a  continuous  solution  applies  in 
the  present  situation,  use  is  made  of  the  fact,  stated  in  the 

t  i 

previous  paragraph,  that  and  s^  have  the  same  sign.  In 
the  ase  of  unloading  from  negative  (compressive)  values 

<  0  SK  <  0  (A0) 

In  addition,  unloading  implies  that  the  derivatives  of  the 
absolute  values  satisfy 

Ji|'  >  0  Kl'  >  0  (A1) 

The  above  inequalities  can  now  be  used  to  see  if  Eq .  (21) 

permits  a  continuous  solution.  Taking  the  derivative  of 
this  equation  with  respect  to  U  one  finds  the  condition 

6U  =  3k'  +  4g'  >  0  (42) 

The  inequality  follows  from  the  fact  that  U  is  inherently 
positive.  From  Eq .  (23) 

iff  <*„»>  ■  -vi  >  0  <43> 

/3  i  i 

Eq.  (25)  gives,  for  T  =  0,  the  relation  gUN,  =  +  —  R4  sx  ~ 

2 

-  -  g6J1  ,  so  that 

dU  (SUN)  =  ^2  84 | S  x |  "  85J1  ~  2R6J1J1  (AA) 
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While  the  first  two  terms  are  positive  the  last  one  is  not. 
Combining  the  last  three  equations  the  continuous  solution 
will  apply  i f 

2/3  g4|sx|'  -  (3k4  +  4g5  +  8g&J1)  >  0  (45) 

The  first  term  is  inherently  positive.  The  second  term  is 
positive  if  the  expression  in  parentheses  is  positive,  but 
as  this  depends  on  the  relative  magnitude  of  8g&J^  <  0 
compared  to  3k^  +  4g^  no  general  statement  can  be  made. 

The  above  relation  investigates  the  situation  on  a  purely 
mathematical  basis.  The  mathematical  condition,  F.q  .  (42), 
however,  expresses  simply  the  fact  that  the  unloading  stress- 
strain  diagram  in  Fig.  2  has  proper  curvature  as  shown.  If, 
therefore,  the  value  of  the  material  constants  give  a  curve 
of  this  type  the  mathematical  condition,  Eq  .  (45),  is  in¬ 
herently  satisfied  and  need  not  be  investigated.  If  the 
stress-strain  diagram  is  of  this  type  it  is  clear,  again  on 
physical  grounds,  chat  no  (unloading)  shock  may  exist  and 
the  continuous  solution  of  Eqs.  (39)  and  (42)  will  be  the 

>V) 

only  solution 

To  obtain  specific  solutions  when  the  general  relations, 

Eqs.  (23)  and  (25),  are  used,  requires  numerical  solution  of 
the  two  Eqs.  (39)  and  (42).  However,  the  principle  may  be 
demonstrated  in  the  special  case  g,.  =  g^  =  0  where  a  closed 
solution  can  be  obtained.  In  this  case  g  is  a  function  of  s 
only,  while  k  is  a  function  of  only  so  that  Eq .  (’9)  can  be 

*7 


Continuous,  but  inappropriate  solutions  of  the  type  dis¬ 
cussed  in  Case  F  may  occur. 
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integrated.  With  b  =  -2/1  g^/9k^ 
g(sx)  =  [Ck(J1)]b 

where  C  is  an  open  constant.  This  relation  must  apply  at 

the  head  of  the  unloading  wave  U  =  U,  where  s  =  s  and 

1  x  xo 

,  which  permits  the  determination  of  the  constant  C 


C 


k(Jl0> 


1/b 


Substitution  of  this  relation  into  the  determinantal  equation 
gives 


3U 


3k ( J  ^ )  +  4g(sxo> 


k(J,) 


lb 


Equation  (46)  may  then  be  usee  to  determine  s 


s  = 
x 


b 

2 

8(sxo> 

‘  k(Jx) 

“  S7 

4 

Si  s4 

L  k(jio) j 

Knowledge  of  the  relation  between  U  and  and  s^.  and 
permits  a  numerical  determination  of  0^  as  a  function  of  U 
in  the  region  to  V  ^  where  varies.  The  typical  shape 
of  this  stress  profile  is  shown  in  Fig.  6a. 


Case  D  -  Pressure  Waves.  Loading. 

Consider  the  case  of  a  stressless  half-space  when  at  the 
time  t  =  0  a  normal  surface  pressure  0  <  0  is  applied,  while 

the  shear  stress  vanishes,  T  =  0.  In  terms  of  the  variable  U 


(46) 


(47) 


(48) 


(49) 
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the  boundary  conditions  are 


0 

x 


(o)  =  a 


xo 


(50) 


Ji(oo)  =  0 

ss(«>)  =  0 


>  (51) 


The  procedure  for  the  determination  of  continuous 
solutions,  if  any,  is  in  principle  the  same  as  in  Case  C, 
except  that  the  expressions  gLD  and  k^  have  to  be  used. 

In  regions  of  continuous  solution  Eqs.  (18),  (19''  and  the 
de te rminan tal  equation 

30 2  -  4g  +  3k  (52) 


must  hold. 

The  possibility  of  continuous  solutions  for  the  present 
case  exists  only  if  the  Eq.  (52)  indicates  that  U  decreases 

1  f 

when  J,  <  0,  s  <0,  while  J,  <  0,  s  <0.  Substitution  of 
l  x  l  x 

Eqs.  (22)  and  (24)  shows  that  the  term  k  in  Eq .  (52)  has 
exactly  the  opposite  behavior,  while  the  second  and  fourth 
terms  of  g  in  Eq.  (24)  change  in  the  direction  required  for 
continuous  solutions.  For  high  stress  levels  the  behavior 
of  k  controls,  but  no  full  statement  can  be  made  without  in¬ 
troducing  specific  numbers. 

As  an  alternative  to  a  computational  investigation  on 
the  basis  of  Eqs.  (22)  and  (24)  qualitative  predictions  can 
be  made  based  on  the  character  of  the  uniaxial  stress-strain 


diagram  implied  in  these  equations.  The  latter  can  be  used 
to  find  a  uniaxial  stress-strain  diagram  which  is  an  approxi¬ 
mation  of  the  actual  diagram,  Fig.  2.  This  approximate 
diagram,  Fig.  5,  ought  to,  and  will  usually  have  the  same 
character  as  Fig.  2,  i.e.,  it  will  show  a  weakening  of  the 
material  at  low  stress  ranges  up  to  point  C,  a  hardening 
thereafter.  The  nature  of  the  pressure  waves  for  this 
situation  is  well  known,  qualitatively.  Numerical  details 
for  specific  values  of  the  coefficients,  F.as.  (22),  (24) 
could  be  obtained  from  Eqs.  (18),  (19)  and  (52). 

a.  If  the  value  of  j°xoj  ■LS  less  than  j°cj  -n  Fig.  5 

the  solution  is  continuous,  changing  between 

points  U ^  and  Uj  ,  .  ig.  6b,  while  the  stress  is 

constant  C  =0  for  U  <  U„  . 

x  xo  —  2 

b.  If  the  value  of  |°xoj>  |  °A  |  ’  is  ■5-arSer  than  J  J  in 
Fig.  5,  the  result  will  be  a  shock,  Fig.  6c.  The 
velocity  is  defined  by  the  slope  of  the  line  OA  in 
Fig.  5.  The  limiting  value  |oTj  is  defined  by  the 
statement  that  OT  in  Fig.  5  is  tangent  to  the  stress- 
strain  diagram  at  0. 

c.  For  intermediate  stress  levels  j°cJ  <  j 0 j  <  |°t| 
there  will  be  a  continuous  precursor,  followed  by 

a  shock.  The  stress  level  up  to  which  the  solution 
Is  continuous  Is  defined  by  point  D  in  Fig.  5,  while 
the  shock  velocity  Is  defined  by  the  slope  of  the 
line  Uli.  iiio  result  has  Che  character  of  Fig.  6d. 


V 


--  17  — 

Case  E  -  Combined  1 ressure  and  Shear. 

The  situation  when  the  applied  surface  load  is  a  com¬ 
bination  of  pressure  and  shear  can  be  discussed  in  a  general 
way,  and  reduced  to  a  combination  of  the  solutions  for 
pressure  and  shear  alone. 

In  the  case  now  considered  where  and  T  must  change, 
all  three  differential  equations  (17),  (18)  and  (19)  must  be 
considered  in  regions  where  the  solution  changes  continuously . 
However,  the  earlier  conclusions  still  hold  that  T  "an  only 
change  if  the  determinant,  Eq .  (20),  vanishes,  while  ,  s^ 
and  thus  can  only  change  if  Eq .  (21)  is  satisfied. 

Potential  solutions  can  therefore  be  found  in  which  changes 
only  say  from  to  >  while  T  changes  only  for  to  , 
the  determination  of  the  distribution  remains  as  outlined  in 
Cases  A  to  D. 

There  is,  however,  also  the  additional  potential  possi¬ 
bility  that  Eqs.  (20)  and  (21) 

U2  =  g  (53) 

U2  =  g  +  k  (54) 

are  satisfied  simultaneously  in  a  region,  to  .  Noting 
that  g  and  k  are  inherently  positive,  regardless  whether  the 
expressions  g  and  k  for  loading,  unloading  or  reloading  are 
used,  the  two  conditions  can  not  be  simultaneously  satisfied. 
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There  can  thus  be  no  region  where  all  three  quantities  , 
and  T  change  simultaneously.  The  pieces  for  const! action 
of  the  solution  are  thus  regions  of  uniform  stress,  regions 
where  T  only  changes,  regions  where  ,  s^_  and  thus  a ^ 
change,  and  compressive  shocks. 


The  sequence  of  the  combination  of  solutions  is  essentially 
the  same  for  all  situations,  and  only  one  case  needs  be  dis¬ 
cussed.  Consider  a  half-space  in  which  at  t  <  0,  the  stresses 

are  J,  =  J,  <  0,  s  =  s  <  0,  T  =  T  everywhere.  At  the 
1  lo  x  xo  *  o  J 

instant  t  =  0  the  surface  loads  are  removed.  The  boundary 
conditions  for  t  >  0  are  then 


O  (0) 

X 

=  0 

T  (0) 

-  0 

J  ,  (°°) 

”  J 

1 

lo 

s  (°o) 

=  s 

X 

xo 

T(w) 

=  T 

O 

Due  to  the  fact  that  in  the  region  of  change  of  t  only 

J2  changes  while  remains  unaffected,  the  construction  of 

a  solution  can  be  started  with  the  region  of  compressive 

change  >  >  ^4  •  This  solution  is  simply  the 

solution  to  the  problem  of  pressure  relief,  when  the  stress 

T  on  the  surface  would  be  maintained,  i.e.,  Case  C  treated 
o 

before.  Let  the  values  of  and  s^  at  point  U ^  found  in 
this  manner  be  J.,  and  s  „  . 


(55) 


(56) 


The  changes  in  stress  which  will  occur  between 


and  <  ^2  ,  are  then  found  for  the  secondary  problem 


Ox('0  =  T(0)  =  0 


(57) 


J,  (“)  =  J 
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1  (to)  =  +T^ 

v2l }  4  Sx2  ^  To 


r(~)  =  t 


(58) 


The  typical  history  of  the  stress  is  shown  in  Fig.  7, 


It  is  not  possible  to  start  the  solution  by  determination 
of  the  region  of  change  in  shear,  and  find  the  pressure  change 
later.  The  reason  is  that  the  solution  must  be  begun  from  a 
boundary  where  all  quantities  J.^  ,  s^  ,  T  are  known.  This 
is  only  true  fcr  U  =  »,  Eqs.  (56).  At  U  =  0  only  the  value  of 
a  =  4  J  +  s  is  known,  not  J1  and  s  separately. 

X  j  1  X  XX 

Other  >  ases  lead  to  appropriate  combinations  of  the 
solutions  shown  in  Figs.  4  and  6,  which  include  pressure 
shocks. 


The  reduction  of  the  problem  of  combined  pressure  and 
shear  to  a  combination  of  the  individual  cases  implies  that 
uniqueness  and  existence  found  in  Cases  A-D,  automatically 
apply  lure. 
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Case  F  -  Elimination  of  Multiple  Solutions. 

In  footnotes  in  Cases  A ,  B  and  C  mention  was  made  of 
additional  "inappropriate"  continuous  solutions.  The  following 
will  give  a  typical  example  and  show  that  such  solutions,  while 
mathematically  correct,  should  be  dismissed  on  physical  grounds. 
This  dismissal  is  based  on  the  requirement  that  any  solution 
containing  a  singularity,  in  this  case  a  step  discontinuity  in 
the  load  at  x=t=0,  must  be  obtainable  by  a  limiting  process 
from  a  rapidly  but  monotonically  changing  surface  load. 

Consider  a  variation  of  Case  B  of  unloading  in  shear, 
where  the  shear  stress  on  the  surface  is  not  completely  removed, 
but  reduced  from  to  niQ  ,  where  0  <  n  <  1.  The  boundary  con¬ 
ditions  are  then 


T(0)  =  nto 
T(“)  =  T 

o 


(59) 


These  boundary  conditions  can  be  satisfied  by  using  Eq .  (34) 

for  unloading  from  t q  to  nTQ  ,  resulting  in  the  monotonic 

stress  history  ABC  shown  in  Fig.  8a.  However,  it  is  also 

possible  to  use  Eq .  (34)  beyond  point  B  reducing  the  shear 

stress  to  a  lower  level,  n,T  <  nt  ,  point  B,  in  Fig.  8b. 

1  O  O  i. 

The  surface  boundary  condition  can  still  be  satisfied  by 
obtaining  from  Eq .  (31)  a  continuous  solution  where  T  rises 
to  the  required  value  nifl  as  indicated  in  Fig.  8b,  points  C^C. 


Points  CC^  for  a  continuous  solution  with  a  rise  in  stress  are 
in  the  required  location,  closer  to  the  surface  than  point  13^ 
because  g  <  applies  for  reloading. 

In  the  alternative  solution  the  time  history  of  the  stress 

at  any  point  is  not  monotonic.  Without  going  into  details,  the 

concept  of  characteristics  requires  that  the  solution  to  the 

problem  of  a  rapid  but  gradual  change  of  the  surface  stress 

from  T  to  nT  would  result  in  a  monotonic  stress  history.  The 
o  o 

solution  shown  in  Pig.  8b  is  therefore  to  be  dismissed.  If 
equivalent  nonmonotonic  results  can  be  constructed  for  other 
problems  they  are  similarly  inappropriate. 


IV  CONCLUSIONS. 


The  investigation  demonstrated  that  in  the  simple  situ¬ 
ations  considered  the  use  of  the  hysteretic  stress-strain 
relations  proposed  in  Ref.  [3]  does  not  lead  to  difficulties, 
i.e.,  solutions  exist  in  all  cases  and  are  unique.  In  some 
cases  uniqueness  is  obtained  only  by  using  an  additional 
physical  consideration  outlined  in  Section  III,  Case  F.  It 
leads  to  the  conclusion  that  solutions  for  the  problem  stated 
here  are  admissable  only  if  the  stress  changes  are  monotonic. 

In  the  case  of  elastic-plastic  materials  a  similar  study 
of  plane  wave  propagation  was  extended  to  the  case  of  the 
superseismic  steady-state  solution  for  a  step  load  progressing 
on  a  half-space.  It  is  expected  that  the  present  work  can 
similarly  be  extended,  thereby  providing  solutions  suitable 
to  check  purely  numerical  approaches. 
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